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THEORETICAL STUDIES OF SOLAR-PUMPED LASERS 


By 

Wynford L. Harries 1 and Zeng - Shevan Fong 2 
SUMMARY 

A method of pumping a CO 2 laser by a hot cavity has been demonstrat- 
ed (ref. 1). The cavity, heated by solar radiation, should increase the 
efficiency of solar pumped lasers used for energy conversion. Kinetic mod- 
eling is used to examine the behavior of such a C0 2 laser. The kinetic 
equations are solved numerically vs. time and, in addition, steady state 
solutions are obtained analytically. The effect of gas heating filling the 

lower laser level is included. The output power and laser efficiency are 

j 

obtained as functions of black body temperature and gas ratios (C0 2 -He-Ar) 
and pressures. The values are compared with experimental results. 

1. INTRODUCTION 

The possibility of using solar collectors on orbiting space vehicles 
and transmitting the energy to earth via laser beams has been discussed 
already (refs. 2, 3, 4). A critical factor for energy conversion is the 
laser efficiency. One suggestion for increasing the efficiency is to heat a 
black body cavity by solar radiation and then let the cavity pump the laser, 
according to the Planck radiation law. If the laser mediun absorbed a nar- 
row bandwidth from the spectrum then the black body cavity would immediately 

refill any "hole" in the distribution. The pumping radiation is thus 

— — - ■ — — - - - 
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replenished and the potential efficiency should be much higher than if the 
laser is pumped directly. 

It is the purpose of this paper to examine the kinetics of a black body 
pumped laser. In addition to directly pumped C0 2 lasers, there are schemes 
whereby N 2 is pumped and the energy transferred to C0 2 , but these are not 
considered here. The direct CC^ lasers usually have He and Ar as additives 
which perform the functions of deexcitation of the lower level and heat 
conduction; these are included. 

The physical mechanisms of the laser are discussed in Section 2. For 
the first time to our knowledge the effect of gas heating, which serves to 
fill the lower laser level and inhibit lasing, is included. The temperature 
rise is calculated, and an estimate of the effect of filling the lower laser 
level from Boltznann statistics is made. 

In Section 3, time varying solutions of the kinetic equations show the 
behavior of the upper and lower laser level densities, and also the photon 

• I 

density and output power vs. time. The solutions show that steady state 
conditions can be achieved. 

Analytical solutions are obtained for the steady state in Section 4. 

The steady state output power is shown as a function of black body temper- 
ature and gas constituent pressures. Estimates of the laser efficiency as a 
function of the variables are also given. 

j. r 

Throughout the discussion it is assumed that the laser tube of 0.4 cm 
radius has a length of 50 cm heated in an oven at 1000-1500 K and contains 
COg+He+Ar at total pressures of 1 to 12 torr in varying ratios. The radi- 
ation from the oven is passed through nitrogen which is transparent to it 
and enters the laser tube. The purpose of the nitrogen, which is flowed, is 
to cool the laser medium. The measured exit temperature of the nitrogen was 
340 K (ref. 1), and this was assumed to be the wall temperature. The gas 
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temperature is slightly higher than this and is discussed in Section 2.2. 

The above nunbers correspond to experimental conditions in reference 1. 

2. PHYSICAL MECHANISMS IN THE LASER 
A flow diagram of the processes occurring during lasing is shown in 
figure 1. 

2.1 Absorption of Radiation 

The black body radiation is absorbed at 4.256y over a bandwidth of 0.1 
y. The radiation flux per unit length entering the laser tube of radius r 
is E(A) dx 2nr when E(X ) is determined by Planck's radiation formula and dA 
is 0.1 y. The fraction absorbed depends on a , the absorption cross sec- 

di 

tion. The latter can be calculated: <x = A 2 A/4ir 2 Av, where A = 424.6 is the 

Q 

Einstein coefficient for the transition (000)^(001), X = 10. 6u, ax = O.ly 
and cr = l^xlO^cm 2 . However, a measured value of the absorption length 

a 

is 1.56xl0- 3 cm at 1 atmosphere (ref. 5), which yields a _ = 2.41xl0 _17 cm. 

a 

I 

Both values were tried. 

Two possibilities exist for the calculation of the fraction absorbed. 

If the absorption is small, then the fraction absorbed per unit volume can 
be expressed approximately as 

Fi = 2 (1 - exp (-(C0 2 ) ° a 2r)/r (1) 

where the radiation is assumed to cross a distance 2r and the amount per 
unit length is proportional to 2irr and per unit volume, 2^rhr z . If the 
pressure increases, however, the absorption will take place at the outer 
edges of the tube. The excited species can, however, diffuse to the center 
in a time t where t * r 2 /0 and D is the diffusion coefficient. Assuming a 
mean free path of 2.9xl0 _3 cm at 1 torr, a total gas pressure of 1 torr 
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(ref. 6) and a temperature of 340°, then t ~ 3 ms, a time comparable with 
other rates. Neglecting diffusion, a second estimate is obtained by 
assuming the radiation is absorbed at the center of the tube between a depth 
x x and x 2 . Then fraction F absorbed per unit volume is 

F 2 = (exp(-(C0 2 )o a x 1 ) - exp(-(C0 2 )a a x 2 ))/(x 2 -x 1 ) (2) 

where here = 0.35, x 2 = 0.45 cm. The effect of the outer layer of gas 
shielding the center of the tube at high pressures is now included, but a 
one dimensional, not cylindrical, geometry is assumed. 

A number of computer runs were made using F 2 , however, the results did 
not agree with the experiment (ref. 1). The computer results showed a sharp 
cutoff of power at a few torr and, accordingly, the coefficient F x was used 
for the following calculations. 

The number of absorption events per unit volume is then 

E(X)dA 

S = — F x =5.76 x 10 19 F!/(exp(3381/T b )-l) (3) 

hv L 

where T b is the black body temperature, which gives the number of C0 2 (000) 
molecules raised to the upper laser level C0 2 (001) per cc per second. 

2.2 Heating of the Medium 

Experimentally it was realized that unless a nitrogen coolant was flow- 
ed past the laser tube, that lasing ceased (ref. 1). The implication was 
that the medium became heated causing the lower level to be filled. A rough 
estimate of the increase in temperature, and its effect on the density of 
the lower level follows. 
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As the efficiency of the laser would probably not be greater than a few 
percent, it is assumed that all the radiation absorbed goes into heating the 
gases. The rate at which heat energy enters a unit volume, dQ/dt, increases 
with the pressure of CO 2 which is included in Fi, 2 . 

dQ/dt = E(x) dx • Fi/4.2 Cals cm- 3 s -1 (4) 

The heat conduction K of a mixture of gases CO 2 , He, Ar depends on the 
gas ratios and the absolute temperature (ref. 3). It is independent of gas 
pressure. 



K 

= 3.27 x 10- 4 S Vt/273 

■ | 

(5) 

where 

S - (- 

+ 0.71 C0 2 / He + 0.32 Ar/He^ 

(6) 

s \l 

+ 4.55 C0 2 /He + 2.67 Ar/H J 


It is assumed that the heat is conducted over an effective area A 
through an effective distance d to the wall at a temperature T . As the 

W 

conduction is proportional to vr. the conduction equation yields that 
the temperature in the center of the gas T g is: 

T g « (Tw + 1.8 x 10 4 • (dQ/dt)/ (S • A/d)f /3 (7) 

The expression is approximate as a one dimensional geometry is assumed. For 

example, with T = 340 K,^ a black body temperature of 1500 K, a total 
w 

pressure of 10 torr (CO 2 20%, He 20%, Ar 60%), tube radius 0.4 cm, then 
approaches 360 K. 
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2.3 Filling of the Lower Laser Level 
The lower laser level which has a density N-j of C0 2 ( 1 00) atoms, will be 
filled according to Boltzmann statistics which depend on T . However com- 
peting against this will be the depletion processes due to collisions. To 
include the temperature effects they must be expressed as a rate of produc- 
tion of N-| vs. time. It is assumed that C0 2 molecules in the tail of the 
energy distribution, which have energies greater than E-j = 0.1649eV, the 
energy of the lower laser level, can collide with the rest of the C0 2 
population with a finite probability of producing C0 2 (100). The density of 
molecules with sufficient energy is (C0 2 ) exp (-E^/kT) and the rate of 
production R 1 from this process is: 


R-j = (C0 2 ) (C0 2 • exp (-E-j/kTg) °v (8) 

' Here k is Boltzmann's constant, a is the cross section for the reaction, v 
the particle velocity. R-j increases as T g increases. 

2.4 Diffusion Losses 

At low pressures losses of particles by diffusion to the tube walls may 
be important. We assume that the diffusion loss time of C0 2 molecules, t d 
obeys (ref. 7): 


x D = (r/2.405) 2 /DY (9) 

where D is the diffusion coefficient for C0 2 in a mixture of C0 2 , He and Ar, 
(Appendix I) and the 2.405 corresponds to losses in an infinitely long cyl- 
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inder. The quantity y is the probability of deexcitation at the wall and is 
somewhere between 0.2 and 1. The diffusion losses apply equally to all 
excited levels of CC^. 

2.5 The Kinetic Equations 

The CO2 laser (Fig. 1) is a four level system. The absorbed photon 
raises the CO^OOO) level to the upper laser level C0 2 (001) the asymmetric 

stretching mode at 0.2912 eV above ground. Lasing at 10. 6y occurs by a 

transition between the (001) and (100) levels, the latter a symmetric 

stretching mode at 0.1649 eV above ground. 

The upper laser level of density N u is depopulated by collisions with 
CCfe, He, Ar (rate constants Kj, K 2 , K 3 ) to the ground level, and by 

spontaneous emission (Einstein coefficient A). Transitions can also occur 
to the lower laser level (100) (K 7 , Kg, K 12 , K 13 , and the spontaneous 
coefficient A^) and the (010) level (K 8 , K 9 , K 14 , K 15 ). 

The (100) level (density N-j ) has a (020) bending mode which is, 

0.0005 eV higher in energy. The conversion (020)^(100) is very rapid, and 
the two are lumped together here. 

The lower laser level is filled by collisions from the upper level (K 7 , 

Kg, K12 » K13) and by the spontaneous emission A u -| . It is also filled during 
lasing by stimulated transitions of rate constant r. The lower level is 

depleted by collisions with C0 2 (K 10 , Kn) and with He (Ki 6 ) > (no data on Ar 
was available), and the transitions end either on the (010) level or the 
(000). However, at high gas temperatures, the lower laser level will be 
filled at a rate R-| . 

The (010) level of density N m is in turn depleted by collision with 
COj, , He, Ar (K 4 , K 5 K 6 ). The reactions are given in more detail in Table 1, 
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with values from four different sources. In the column under reference 9, 
the values are calculated from the expression 

- 1/3 - 2/3 

In K = A + BT + CT (10) 


where the constants A, B, C, are given for the different reactions. The 
values may not be accurate when extrapolated to 360K. References 8, 10 and 
11 give experimental values. 

The unknown quantities are then assumed to be N , N-j , N m and N the 
10. 6y photon density. In the following equations, the quantities in paren- 
thesis (CO^) etc. are the densities of the species. The rate equations for 
the four species are: 



dt 



[(*! + K 7 + Kg + K 9 )(C0 2 ) + (K 2 + K 12 + K 14 )(He) 


+ (K 3 + K 13 + K 15 )Ar] - N u (A + A ul )- r - N u/ (11) 


dN 1 

— - r + N u A ul + N u [(K 7 + K 9 )(C0 2 ) + K 12 (He) + Ki 3 (Ar)] 
dt 

- [(Kiq + K 11 )(C0 2 ) + K 16 (He) ] - N e/ + R ] (12) 


dN m 

— = N 1 (K 10 + 2K 11 )(C0 2 ) + N u [(K 8 +K 9 )(C0 2 ) + K 14 (He)+K 15 (Ar)] 

- N m [K 4 (C0 2 ) + K 5 (He) + K 6 (Ar)] - N m/ (13) 

t D 
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dN 


(14) 


dt 


+ G 


A ul N u 


7t 


The stimulated emission rate r is given by: 


T = N a e c [N u - N-j ] (15) 

and is proportional to the stimulated emission cross section whose value 

4 

was calculated: <x - \ n A ,/ 4tt 2 c AA , where A = 10. 6u, A . = 0.186 and 

e e ul e e ul 

AA g was taken as the bandwidth from Doppler Broadening. Here c is the 
velocity of light. Assuming a gas temperature of 360°, then AA g = 0.217A 
and o = 9.14xl0“ 17 cm 2 . 

£ ' I 

i ■ 

In equation 15 the quantity G is a geometrical factor which takes 

into account that the spontaneous emission is isotropic; G = 2r 2 /L 2 where 

r and L are the radius and length of the laser. The quantity x c is the 

containment time for photons, t = -2L/c In (r x r 2 ) where r x = 0.99, r 2 = 

c 

0.95 are the reflectivities of the laser mirrors. 

The equations can be simplified by noting that N^, <<(C0 2 ), and 

that (CO 2 ) can be regarded as constant. Then the following constants can be 
defined: i 


ai = [ (Ki +K 7 +K 8 +K 9 ) (CO 2 ) + (K 2 + K 12 + Km) (He) + (K 3 + K n + Kis)(Ar) 


+ A + A . + 
ul 



] 
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1 

a 2 = [ (Kio+Kn) (C0 2 ) + Ki 6 (He) + — )] 

t d 

A' = [ {K 7 +K 9 ) (C0 2 ) + K i 2 (He) + K 13 { Ar) + A ul ] 
a 3 = [K 4 (C0 2 ) + K 5 (He) + K 6 (Ar)j 
a* = [(K10+ 2 Ki 1 ) (C0 2 ) ] 


a 5 = [(K 8 +K 9 )(C0 2 ) + K 14 (He) + Ki S (Ar)] 


Then equations 11 through 14 can be reduced to 


dN. 


dt 


S - aiN. 


dN, 


dt 


+ A'N U - a 2 N 1 + R ] 


dN 


m 


— = asN u + ^ N 1 ' a 3 N m 


dN 


dt 


= r + SA ol N u - n /t c 


(16) 


(17) 


(18) 


(19) 
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3. TIME VARYING SOLUTIONS 

Computer graphic solutions of equations 15 through 19 are shown in 
Figure 2. It was assumed for convenience that the pumping signal from the 
black body rose from zero to full value in a time of 100 us, and that the 
signal was proportional to sin 2 2n (t/10“ 4 ) (where t is the time in sec- 
onds) up to 100 us; thereafter it was constant. The pumping signal is shown 
in the upper curve. 

In the second plot, N u , N-j and (N u -N-|)xlO are plotted. The densities 
of the upper and lower laser levels increase with time, but at first, when 
the signal is small, N-j > N y . (This can.be seen on amplified plots). How- 
ever, N^ overtakes N^ at 28us, and the signal (N u -N-| ) starts to grow. It 
increases with time until threshold at t = 53 us. Then N u falls and N-| 
rises, and the inversion (N u -N^) falls. At the same time the photon density 
N increases sharply and then falls to zero as the inversion falls. The 
upper level is then increased again until a second light pulse appears and 
so on. This behavior is characteristic of kinetic equations similar to 
these, where is small; if a g is much larger these violet oscillations 
are smoothed out. 

The inversion quickly assumes a constant value, and at 135 us the value 
of the photon density N also becomes constant. This implies that continu- 
ous lasing is possible. 

4. STEADY STATE SOLUTIONS 

As they stand, equations 16 through 19 are nonlinear. In the steady 
state, however, the stimulated emission rate far exceeds the spontaneous 
rate, and in equation 19, r » GA^i N u or: 
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iody pumped laser. The pumping signal 
to sin 2 (irt/2xl0 -4 ) up to lOOys, then 
ddle diagram shows Nu> N] and (Nu-N] )xlO 
the 10. 6u photon. 



r = 



( 20 ) 


and from 15 and 19 


O CT 

e c 


= const 


( 21 ) 


The constancy of the inversion N^-N^ is evident in Figure 2, and the 

value is given by equation 21. Replacing r in equations 16 through 19 by 

N/t the solutions are 
c 


S + R 1 + a 2 / T c * e c 
a^ + &2 — A 1 


n i ■ %- 


TOC 

c e 


( 22 ) 


(23) 


a 5N u + 


N m = 
m 


(24) 


a l 

1 _ 


a l a 2 


■ T c R l a i ' 

dLi+a 2 -k' 


o e c(a 1 +a 2 -A’ ) 


a 1 +a 2 -A‘ 


(25) 


The effect of gas temperature is included in and as R^ increases the 
value of N falls. 


The value of N from equation 25 agrees with the value, late in time, 
from the time varying computer solutions. 
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4.1 Steady State Power Output 

The evaluation of N enables the power output of the laser to be 
evaluated: 

1 

P = — Nc hv T irr 2 (26) 

2 r 

The factor 1/2 takes into account that half the photons are travelling away 
from the output mirror of transmission T^ = (l-r 2 ), c is the velocity of 
light, hv the energy of the lasing photons of 1.868xlO- 20 J. Experiment- 
ally it was found that the output beam had a radius of 0.3 cm corresponding 
to a Gaussian profile in a 4 cm radius tube, (ref. 13) and for T r = 0.05, 

• I 

P = 3.97x10- 12 N (W) (27) 

Plots of power output vs. total pressure p (torr) for different black 
body temperatures are given in Figure 3. The gases were in the ratio C0 2 
1656, He 4%, Ar 80°. The output power rises with black body temperature, but 
as p increases, the power falls. The effect is due to the term in the 
equations; on making R-j = 0 the output power rose with increasing pressure. 

The variation of output power with black body temperature at a fixed 
total pressure of 7 torr, but with the same gas ratios, is shown in Figure 
4. A linear rise is observed. 

The gas temperature of T^ can be altered by assuming different wall 

temperatures T . Assuming a gas mixture as before of 7 torr (16% C0 2 , 4% 
w 

He, 80° Ar) and a black body temperature of 1450K, then the power output 
vs. Tg is plotted in Figure 5. It shows that an increase in gas 
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0 12.5 . 

Total Pressure p (torr) 


Figure 3. Power output of black body pumped laser vs total pressure 
for various black body temperatures. The gases were 
C0 2 -16%, He-4%, Ar-80°. The wall temperature was 340K. 
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Black Body Temperature (K) 


Figure 4: Power output vs black body temperature. Total pressure = 7 torr, 

16% C0 2 , 4% He, 80% Ar, wall temperature 340K. 
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-»■ p (torr) 


Figure 6. Efficiency (as defined in equation 29) of the black body 
pumped laser vs total pressure p for various black body 
temperatures. The gases are C0 2 16%, He 4% and Ar 80°. 
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temperature greatly reduces the output power indicating that cooling of the 
medium is essential. 


4.2 Efficiency 

The input power P^ is the power absorbed by the laser: 

P in « E(X) dX 2nr L (l-exp(-(C0 2 )o a 2r) (28) 

and the efficiency n is defined here as output power/the power absorbed 
from the black body radiation: 


n = P/P. (29) 

in 

I 

(It should be noted that the black body pumped laser has other "efficien- 
cies" besides n, however, these are not considered here.) 

Using equation 29 to define the efficiency. Figure 6 shows plots of n 
vs. p for different black body temperatures. There seem to be high effi- 
ciencies at very low pressures where the output power is small. In the 
regime around 7 torr where experimentally the highest output powers were 
observed, the efficiency is only a few percent. 

CONCLUSIONS 

The results here depend critically on the values of the rate coeffi- 
cients in Table 1. They are insensitive to K 4 , K 5 , K 6 , but are very de- 
pendent on Kg, K n and k 13 • The most sensitive parameter seems to be K n 
where a change of a few percent altered the output power by several times. 
Also t p was a sensitive parameter. 

Unfortunately, the rate coefficient for R-| was not known. The values 
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chosen indicate that about one collision in 10 5 of two ground state atoms 
resulted in the creation of 'one C 02 ( 100 ) which does not seem unreasonable. 

Figures 3, 4, and 5 show reasonable qualitative agreement with 
experimental measurements (ref. 1). However, it did not seem possible to 
get agreement with experimental measurements of power output vs. total 
pressure (ref. 12), at a black body temperature of 1450K where the ratios 
CC^rHerAr were varied. In view of this, these results must be treated with 
caution, and must be regarded as preliminary. 
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APPENDIX 


Diffusion of CO 2 When Mixed with He and Ar 


The mean free paths of molecules in their own gas at 1 torr for CO 2 , He 
and Ar are 2. 7xl0- 3 cm, 13.1xl0 -3 cm and 4.6xl0- 3 cm, respectively (ref. 7). 
These values give the cross sectional diameters, and it is then assumed that 
for a collision between unlike molecules that the cross section diameter is 
the arithmetic mean of the individual diameters. Adding the collision 
frequencies of CO 2 -CO 2 , C0 2 -He, C0 2 -Ar, the total collison frequency v for 
C0 2 is: 


v 


/ He 

2.37xl0 -11 (C0 2 ) ( 1 + 1.79 + 

' C0 2 



(Al) 


The diffusion coefficient D is then 


1 v r 

D = — = 1. 928x1 0 6 T/v 

3 v 


where v = 2405 vr is the random velocity of CO 2 . 
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